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Abstract 

A  scheme  for  ultrahigh-speed  all-optical  switching  based  on  semiconductor 
heterostructures  was  explored  theoretically.  The  approach  was  anticipated 
to  be  applicable  to  systems  operating  in  the  10  to  100  gigabit/sec  regime. 

Our  role  was  to  simulate  the  carrier  dynamics  in  semiconductor  heterostruc¬ 
tures  under  ultrafast  (picosecond  and  subpicosecond)  laser  excitation.  Pulse 
characteristics  were  forund  to  coherently  populate  and  depopulate  the  semi¬ 
conductor  heterostructure  of  carriers  on  the  picosecond  timescale. 
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The  transmission  of  data  through  optical-fiber  communications  systems  is  ultimately  limited 
by  the  optical  bandwidth  of  silica  fiber  well  in  excess  of  1  terahertz.  Current  systems,  however, 
are  limited  to  a  small  faction  of  that  bandwidth  due  to  limitations  associated  with  modulating 
and  demodulating  optical  data  streams.  One  approach  to  high  bandwidth  optical  communications 
systems  is  optical  time-division  multiplexing  (OTDM)  several  data  channels.  Each  individual 
channel  is  modulated  at  rates  consistent  with  high-frequency  electronics.  Putting  the  channels 
together  (multiplexing)  and  then  teasing  them  apart  (demultiplexing)  requires  optical  switches  that 
are  as  fast  as  the  aggregate  bit  rate.  Although  nonlinear  fiber  loop  devices  have  been  demonstrated 
with  the  requisite  speed,  a  purley  semiconductor-based  approach  is  desired  to  facilitate  integration 
of  optical  switching  devices  with  other  components. 

We  explored  the  use  of  specially  shaped  optical  pulses  to  devise  a  possible  switching  scheme 
that  could  operate  at  the  requisite  rates  in  a  purely  semiconductor-based  structure.  Namely,  by 
coherently  driving  carrier  populations  up  and  then  back  down,  we  planned  to  demonstrate  theo¬ 
retically  a  route  to  optical  switching  that  would  meet  the  requisite  100  gigabit /sec  rates,  but  for 
pipeline  architectures.  The  quintessential  phenomenon,  well  known  from  atomic  physics,  is  Rabi 
flopping.  The  main  question  is:  to  what  extent  can  Rabi  flopping  associated  with  the  creation 
and  coherent  destruction  of  excitons  be  controlled?  Whereas  in  atoms  the  approach  is  based  on 
the  optical  Bloch  equations  (which  describe  the  coherent  coupling  of  a  two-level  system  to  an  elec¬ 
tromagnetic  field),  for  semiconductors  the  appropriate  description  is  based  on  the  semiconductor 
Bloch  equations  which  describe  the  evolution  of  the  interband  polarization  and  the  electron  and 
hole  populations  coupled  to  the  electromagnetic  field. 

Our  most  important  result  is  that  one  can  use  specially  temporally  shaped  picosecond  optical 
pulses  to  coherently  control  carrier  populations  in  quantum  wells.  In  particular,  strongly  chirped 
and  frequency-modulated  pulses  were  found  to  lead  to  pronounced  population  trapping  effects  and 
excitonic  Rabi  oscillations.  Despite  the  fact  that  these  effects  were  strong,  it  is  not  clear  that  they 
are  very  useful  directly  for  optical  switching,  since  following  trapping,  the  carrier  population  goes 
into  free-carrier  states  (i.e.,  not  back  down  to  the  crystal  ground  state). 

Due  to  the  termination  of  this  program  well  before  the  goals  could  be  reasonably  met,  portions 
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of  this  work  were  carried  out  under  other  funding  subsequent  to  termination.  Nevertheless,  it 
should  be  noted  that  the  funding  from  AFOSR  was  essential  to  begin  work  in  this  area  in  three 
key  ways.  (1)  It  funded  Dr.  Harshawardhan  Wanare  as  a  post-doctoral  research  fellow  for  less  than 
a  year.  Dr.  Wanare  is  an  expert  at  coherent  atom/field  interactions;  we  required  his  expertise  to 
carry  over  the  relevant  physics  to  the  semiconductor  domain.  (2)  Due  to  the  widespread  and  well 
recognized  difficulties  finding  post-doctoral  research  fellows  in  optics,  it  was  not  until  April  1998 
that  Dr.  Wanare  could  be  hired.  Nevertheless,  it  became  evident  that  the  numerical  problems  we 
were  about  to  tackle  were  considerably  more  difficult  than  we  had  anticipated.  We  therefore  used 
funds  from  this  grant  to  purchase  computer  equipment  (see  attached  memos).  This  equipment  was 
required  to  perform  computations  supported  by  this  grant  [D.  S.  Citrin  and  W.  Harshawardhan, 
“Terahertz  sideband  generation  in  quantum  wells  viewed  as  resonant  photon  tunneling  through 
a  time-dependent  barrier:  An  exactly  solvable  model,”  Phys.  Rev.  B  60,  1759  (1999);  W.  Har¬ 
shawardhan,  S.  Hughes,  and  D.  S.  Citrin,  “Ultrafast  coherence  effects  in  semiconductor  quantum 
wells  using  modulated  optical  fields,”  Centenniel  Meeting  of  the  American  Physical  Society,  Talk 
IC15.07]  as  well  as  to  carry  out  work  initially  anticipated  to  be  funded,  but  ultimately  not  funded  by 
this  grant  but  supported  by  other  sources  as  an  emergency  measure  [S.  Hughes,  W.  Harshawardhan, 
and  D.  S.  Citrin,  “Excitonic  state  trapping  and  quasiadiabatic  population  transfer  in  a  two-band 
semiconductor,”  Phys.  Rev.  B  60,  15523  (1999)].  (3)  The  grant  partially  supported  two  Ph.D. 
students  (Michael  Feise  and  Alex  Maslov)  who  contributed  to  the  numerical  computations. 

In  addition,  this  work  has  had  a  spin-off  which  we  have  pursued  subsequently  under  different 
funding.  In  S.  Hughes  and  D.  S.  Citrin,  “Extremely  tunable  terahertz  emission:  Coherent  popula¬ 
tion  flopping  in  a  dc-biased  quantum  well,”  Optics  in  1999 ,  Optics  and  Photonics  News  10  (12), 
44  (1999).  In  this  work  we  showed  theoretically  that  it  may  be  possible  to  exploit  excitonic  Rabi 
flopping  in  semiconductor  quantum  wells  as  a  source  of  coherent  transient  electromagnetic  radia¬ 
tion  in  the  1-10  THz  regime.  This  is  exciting  since  there  are  few  convenient  sources  of  coherent 
transient  radiation  in  this  portion  of  the  spectrum. 
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DEPARTMENT  OF  PHYSICS 
Washington  State  University 


AFOSR,  2447-0020 
Grant  F496209710535 


May  12,  1999 


PULLMAN,  WA  99164-2814 

TELEPHONE:  (509)  335-3698 
TELEFAX:  (509)  335-7816 
citrin@wsu.edu 


To  Whom  It  May  Concern: 

This  letter  is  to  justify  the  expenditure  of  45  %  of  grant  F496209710535  for  equip¬ 
ment  when  no  equipment  was  included  in  the  original  budget. 

The  start  date  of  this  three-year  project  was  11/1/97.  The  budget  included  the 
salary  for  a  postdoctoral  research  associate  for  the  full  three  years.  However,  it  was 
not  until  4/27/98  that  I  was  able  to  hire  a  post-doctoral  research  associate  to  work 
on  the  project.  When  it  had  become  evident  that  there  would  be  a  delay  in  hiring 
the  post-doctoral  research  associate  and  that  the  computational  complexity  of  the 
project  exceeded  our  expectations,  I  decided  to  use  some  of  the  funds  that  were 
originally  budgeted  for  salaries  (plus  associated  fringe  benefits  and  indirect  costs)  to 
purchase  computer  equipment.  I  believed  that  this  was  essential  to  meet  the  long¬ 
term  goals  of  the  project.  This  included  three  Digital  workstations  (DEC  DPW-500) 
as  well  as  associated  peripherals.  The  original  requisition  for  the  equipment  dates 
from  11  December  1997  whereas  the  purchase  order  is  dated  26  February  1998. 
The  equipment  itself  was  received  during  June  of  1998,  but  was  not  paid  for  until 
much  later  because  the  vendor  did  not  fulfill  all  the  terms  of  the  agreement  until 
then.  These  machines  were  used  by  me,  my  graduate  students,  and  the  postdoctoral 
research  associate  (after  4/27/98)  to  carry  out  nonlinear  optical-pulse  propagation 
in  low-dimensional  semiconductor  structures  (microcavities),  as  had  been  outlined  in 
my  proposal.  A  technical  report  was  submitted  detailing  the  progress  made  during 
the  first  year. 

However,  on  August  20,  1999,  I  received  notification  from  Dr.  Alan  Craig  that 
the  project  would  not  be  continued  past  the  first  year,  which  ended  10/31/98  (see 
copy  of  the  notification  attached).  This  happened  because  the  Innovative  Science 
and  Technology  directorate  at  BMDO  was  "disestablished”  and  was  not  due  to  any 
failure  on  my  part.  In  spite  of  the  fact  that  it  is  currently  unfunded,  I  have  continued 
to  work  on  the  project  using  the  computer  equipment  which  was  purchased  by  the 
project.  A  list  of  publications  Submitted  and  conference  presentations  given  that 
were  supported  by  this  project  is  available. 

Please  do  not  hesitate  to  contact  me  for  further  information. 


David  S.  Citrin 

Assistant  Professor  of  Physics  and  Materials  Science 


From:  alan.craig@afosr.af.mil 
Date:  Thu  Aug  20  15:20:30  1998 
To:  <citrin@wsu.edu> 

Subject:  BMDO  funds  eliminated 

Dear  David, 

The  Innovative  Science  and  Technology  directorate  at  BMDO  has  been 
disestablished.  (Lou  Lome,  the  IS&T  champion  of  efforts  to  explore  high 
bandwidth  optical  networks  for  distributed  computation  with  whom  I  have 
worked  for  many  years,  has  been  reassigned  to  a  Simulation  and  Modeling 
group.)  As  a  result,  funds  previously  allocated  to  agents  acting  for  BMDO 
to  pursue  technologies  that  support  its  goals  have  disappeared  and  programs 
dependent  on  them  have  been  discontinued.  Unfortunately,  the  termination  of 
the  IS&T  function  affects  the  continuation  of  your  grant  through  my  office. 

I  have  struggled  to  ameliorate  the  impact  of  this  disaster,  but  am  unable 
to  discover  any  means  to  circumvent  the  budgetary  difficulties.  As  a  result, 
the  remaining  option  increments  on  your  grant  will  not  be  exercised. 

Your  investigation  into  the  theory  of  coherent  electronic  states  in 
semiconductors  has  only  just  begun  in  this  program,  so  I  apologize  for  the 
brief  support.  Together  we  gambled  that  the  BMDO  IS&T  program  would 
survive;  now  we  endure  the  awkwardness  of  losing  the  bet.  Only  because  I 
hope  so  fervently  that  perceptions  in  this  arena  will  yield  new  technology 
did  the  bet  make  sense;  now,  in  a  program  that  supports  primarily 
experimental  work,  it  can  no  longer  be  justified.  With  time  and  persistence 
this  area  should  garner  renewed  support.  Perhaps  your  efforts  over  the  past 
year  have  sufficiently  established  this  field  that  it  will  attract 
recognition  and  investment  at  agencies  with  more  expansive  vision. 

Thank  you  for  your  participation  in  the  AF0SR/BMD0  program. 

Sincerely, 

Alan 

August  13,  1998 


Alan  E .  Craig 
AFOSR/NE 

110  Duncan  Avenue,  Suite  B-115 
Bolling  AFB,  DC  20332-8050 
Tel:  (202) -767-4934 

Fax:  (202) -767-4986 

E-mail:  alan.craig@afosr.af.mil 
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The  propagation  of  light  through  quantum  wells  in  the  vicinity  of  an  excitonic  resonance  subjected  to  a  THz 
electric  field  is  considered.  The  treatment  makes  explicit  the  relationship  between  THz  sideband  generation  and 
resonant  photon  tunneling  through  a  time-dependent  resonant-tunneling  structure.  A  method  for  the  coherent 
control  of  excitons  in  quantum  wells  is  also  suggested.  [SOI  63-1 829(99)1 1927-1] 


The  basis  for  much  work  on  transport  through  harmoni¬ 
cally  varying  time-dependent  potentials  dates  back  to  the 
work  of  Tien  and  Gordon1  on  tunneling  between  supercon¬ 
ducting  films.  Subsequent  developments  in  semiconductor- 
heterostructure  growth  and  split-gate  devices  largely  account 
for  the  renaissance  of  interest  in  this  area.  A  parallel  devel¬ 
opment  is  the  ffee-electron  laser.  The  modulation  via  micro- 
wave  and  THz  fields  of  transport  properties  of  quantum- 
point  contacts,2  single-electron  devices,3  and  resonant¬ 
tunneling  structures4,*has  been  studied.  Effort  has  also  been 
put  into  observing  THz  sidebands  (TS’s)  in  the  optical  spec¬ 
tra  of  semiconductor  quantum  wells  (QW’s).6  Despite  sub¬ 
stantial  theoretical  progress  in  the  area  of  tunneling  through 
time-dependent  potentials,  experimental  results  have  been 
scant.  And  apart  from  the  common  role  of  the  low-frequency 
field  (microwave,  THz)  in  inducing  mutiphoton  transitions  in 
the  system,7  an  explicit  elucidation  of  the  connection  be¬ 
tween  resonant  tunneling  through  time-dependent  potentials 
and  TS  generation  (TSG)  was  hitherto  unavailable.8  The  con¬ 
nection  though  may  prove  valuable  since  resonant-tunneling 
devices  based  on  semiconductor  heterostructures  introduce 
numerous  complications  prohibiting  the  controlled  variation 
of  a  single  parameter.  These  difficulties  have  hampered  a 
detailed  analysis  of  transport  experiments  through  time- 
dependent  potentials.  Typically,  experimental  results  have 
been  modeled  using  Tien-Gordon  theory  which  only  captures 
the  most  basic  effects.7,9  Recently,  however,  exciting  work 
has  been  published  showing  that  the  quantum  mechanical 
phase  of  electrons  can  be  measured  in  tour  deforce  transport 
experiments10,11  lending  urgency  to  the  exploitation  of  opti¬ 
cal  analogues  where  phase  measurements  are  routine.  Optical 


analogues  of  time-dependent  transport  problems  therefore 
hold  out  the  intriguing  prospect  of  realizing  a  number  of 
effects  that  have  been  theoretically  predicted.  Already,  con¬ 
siderable  experimental  work  on  superluminality  and  its  im¬ 
pact  on  issues  concerning  tunneling  times  has  been  carried 
out.12 

In  this  study  we  present  a  theoretical  treatment  of  light 
propagation  through  a  QW  in  which  the  exciton  energy,  line- 
width,  and  oscillator  strength  are  modulated  with  bandwidths 
in  the  100-GHz  to  5-THz  range  (henceforth  denoted  THz 
field),  such  as  provided  by  excitonic  Stark  effects  induced  by 
a  THz  field  F (t)  polarized  in  the  QW  plane  x  or  the 
quantum-confined  Stark  effect  with  F (f)  parallel  to  the 
growth  direction  z.  The  formalism  employed  makes  explicit 
the  close  relationship  between  TSG  and  resonant  tunneling 
of  photons  through  the  QW  in  the  presence  of  a  time- 
dependent  dielectric  constant.  As  is  well  known,13  the  wave 
and  Schrodinger  equations  can  be  mapped  from  one  to  the 
other  via  (a>/c)n(<u,r)^{(2w/^2)[£(co)  — V(r)]}1/2,  re¬ 
spectively,  where  the  symbols  here  have  their  usual  mean¬ 
ings.  Previous  studies  of  TSG  (Ref.  6)  have  avoided  the  de¬ 
tails  of  the  propagation  effects  which  bring  to  the  fore  the 
close  analogy  between  the  two  classes  of  phenomena. 

We  consider  linear  optical  propagation  although  the  THz 
field  may  be  strong.  Thus,  we  are  concerned  with  the  low- 
density  regime  in  which  excitation-induced  dephasing  and 
saturation  may  be  neglected.  Let  an  optical  pulse  [amplitude 
/(f)]  be  incident  normally  on  the  QW  in  the  vicinity  of  an 
excitonic  resonance.  The  reflected  R(t )  and  transmitted  T(t) 
optical  fields  are  given  by14 


no' 

0 


-  f_Jt'  x(m') 

PRB  60 


t 


'1(f) 

R(f) 


(1) 


0163-1 829/99/60(3)/!  759(5)/$  1 5 .00 


1759 


©1999  The  American  Physical  Society 


1760 


D.  S.  CITRIN  AND  W.  HARSHAWARDHAN 


PRB  60 


with  *(f,f')  the  optical  susceptibility  (to  be  discussed  be¬ 
low)  and  the  integrals  operate  on  everything  to  the  right.14 
We  eliminate  R(t)  from  Eq.  (1)  to  obtain 

T(t)  =  I(t)  +  (2) 

If  F=0,  then  x( =  and  Eq.  (2)  is  easily  solved 

in  the  frequency  domain;  however,  for  F^O,  the  susceptibil¬ 
ity  no  longer  depends  simply  on  t— 1\  but  explicitly  on  the 
two  time  variables.  Nevertheless,  we  shall  see  that  even  in 
this  case,  under  certain  conditions  an  explicit  solution  ob¬ 
tains. 

The  optical  susceptibility  *(t,r')  in  the  presence  of  F(t) 
can  be  treated  as  in  Ref.  15  or  solved  for  numerically;16 
however,  as  we  show,  a  simple  model  can  be  considered 
which  leads  to  a  closed-form  solution  of  Eq.  (2).  To  motivate 
a  physically  reasonable  ansatz  for  first  consider  the 

case  F=0.  For  a  single  resonance  (e.g.,  retaining  the  Is  ex- 
citon,  but  neglecting  the  other  bound  states  and  continuum) 
we  have14  *(f-f')  =  re~'(€",r)(r_r,)  where  T  is  the  exciton 
radiative  width,17  e  is  the  frequency  of  the  exciton  measured 
from  the  crystal  ground  state,  and  y  is  the  nonradiative  line- 
width.  The  parameters  e,  y,  and  T  are  for  F(r)  =  0.  All  other 
interband  resonances  of  the  QW  are  neglected  and  the  line 
shape  is  assumed  for  convenience  to  be  Lorentzian.  If  the 
Hamiltonian  governing  the  interband  electronic  excitations 
of  the  QW  depends  on  time  through  a  THz-frequency  elec¬ 
tric  field,  but  with  frequency  content  sufficiently  below  any 
intraband  resonance,  we  can  account  for  the  effects  of  the 
THz  field  through  an  adiabatic  Stark  shift,  line  broadening, 
and  modification  of  the  oscillator  strength.  (For  the  param¬ 
eters  investigated  here,  ultrahigh-quality  ZnSe  or  GaN  QW’s 
might  be  the  best  candidates  due  to  their  high  exciton  bind¬ 
ing  energies.)  We  can  thus  use  the  ansatz  x(t>t') 
=  G(t)G*(t')  with  =  where 

e(/),  T (0,  and  y(t)  all  depend  on  time  due  to  the  influence 
of  F.18  The  transfer-matrix  method  is  a  formulation  of  the 
scattering  approach  to  such  problems  and  is  well  known  to 
follow  from  the  tunneling  Hamiltonian  for  the  appropriate 
structure.18  The  reflectivity  and  transmission  coeffeicients 
are  equivalent  to  the  Green  function  of  the  appropriate  tun¬ 
neling  Hamiltonian. 

Equation  (2)  can  be  converted  into  a  differential  equation. 
The  solution  is 

T(t)=I(t)+  P  dt'  (3) 

J  -00 
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FIG.  1.  Reflectivity  |Rw(fc>  +  /x£)|2  int°  harmonic  (jl  normalized 
to  peak  reflectivity  of  the  unmodulated  fundamental  (I7T')2 
=  1.7%  as  a  function  of  detuning  S  from  exciton  line  center.  The 
parameters  are  e1  =  2.5  (solid)  and  1.52  ps-1  (dashed);  in  all  cases 
y0=  1.52  ps-1,  T  =  0.1  ps-1,  and  y!  =  0. 

^o-rjcos^f  where  £=H  if  Fdc>Fac  or  £=2 Cl  if  Fdc 
=  0.  Typical  values  may  be  found  in  Refs.  19.  We  set  e, 
=  £,■— iyt  and  ej  =e,— iT,-  with  (  =  0,1.  If  Fdc=0,  Stark  ef¬ 
fects  are  quadratic  and  the  parameters  are  modulated  at  fre¬ 
quency  2 Cl.  This  leads  to  even  TS’s  of  Cl  as  is  expected  for 
a  system  with  inversion  symmetry.  If  Fdc±  0,  the  system 
lacks  inversion  symmetry  and  both  even  and  odd  TS’s  occur. 
Henceforth,  we  assume  1^  =  0.  Let  l{t)  —  e~l<i)t.  Then 

nt)  =  T„(t)  =  e-iui 

+  r  e~iwt  J  dtf  £f-sin  £*'] 

(4) 

Set  t=  t  — 1\  7=(f  +  r72),  and  sin  £f- sin  £7 

=2  sin(£7/2)cos  £F.  Fourier  transform  Tw(  o) ' ) 

=  S\wlidtelm',TJLt)  to  obtain  r.(«')  =  (2ir/f)[<<» . 

^ot-~ a)' ,fi£ ]  with 


This  defines  the  response  function  S(t—t’)  +  x(td')  with 
X(td')  =  g(t)g*(t'),  and 

g^  =  rm^t)e-U'di'  {e(t')-iW)+ni'))) 

Using  T=I+R  [cf.  Eq.  (1)],  we  have  R(t) 

=  f'-„dt'  xitd'W). 

In  the  following  we  assume  the  THz  field  F(f)  =  Fdc 
+  Fac  cos  fir  is  applied  to  the  QW  with  Fdc  (which  may  be 
zero)  a  dc  bias  and  Fdc||Fac.  We  write  the  modulated  param¬ 
eters  as  e(f)  =  e0  +  £]  cos£f,  y(t)  =  yQ+  yx  cos  £/,  and  Tfr) 


^„(«)  =  2irAS 


1 


*-i  A2-(kC/2)2 


XJ, 


/e  i*\  (st  k 


(k+n)n\^ j 2A"  ’  ^ 


a}  =  (o)  +  o),)f 2  the  average  of  the  incoming  and  outgoing 
frequencies,  /x  an  integer,  A  =  A^(a))  =  a>- e(0=(o- fji£/2 
—  Ef0  the  complex  detuning  between  the  average  frequency 
and  the  radiatively  renormalized  excitonic  resonance,  and 
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FIG.  2.  Total  reflected  probability  /?tot(e)  for  the  same  parameters  as  Fig.  1  (solid).  Note  that  for  clarity,  the  region  |£|>20  ps  1  is 
multiplied  by  a  factor  of  50. 


2^1 1  the  sum  restricted  to  integers  k  of  the  same  parity  as  fx. 
Rw( (o')  is  obtained  from  the  previous  formula  by  replacing 

( (o-E0)l(<o-ef0)SU)  Qjf  by  iT/(a)-eo)^«,«#  •  ^  first  term 

leads  to  the  unmodulated  result,14  e\  ,yj  =  0.  The  other  term 
gives  the  TS’s  as  well  as  the  modified  transmission  of  the 

fundamental.  Finally,  for  e  j  =  e  i  (  y  i  =  0) ,  direct  integration 
gives  for  the  total  reflected  probability  per  unit  time 


RJtu)  =  (C/2ir)fl*,ldt\Rm(t)\2-(C/2v) 
X(r/r')Re/?w(fi)>  (optical  theorem)  where  r'  =  r  +  y0. 

The  foregoing  expressions  for  R  and  T  are  compact;  how¬ 
ever,  it  is  useful  to  reexpress  the  Kapteyn  series  for  JC^o)) 

in  Eq.  (5)  as  a  power  series  in  z  =  /£  (Ref.  20)  to  identify 

the  various  multi-THz-photon  processes  that  contribute  to  a 
given  TS: 


r 


fjb  even 


jjl  odd. 


(6) 


Series  (5)  and  (6)  converge  in  \ze ^  **/(!  +  Jl-z2')\<l 
and  for  arbitrary  2A/£  not  an  even  (odd)  integer.  0  The  terms 
to  a  given  order  in  £j  summed  over  fx  are  the  nonlinear 
susceptibilities.  For  fixed  /x,  the  term  of  degree  p  in  e\  de¬ 
scribes  pth-order  multi-THz-photon  processes  contributing 
to  the  TS.  This  in  turn  allows  one  to  write  the  quantum 
mechanical  transmission  amplitude  as  a  sum  of  nonlinear 
susceptibilities,  as  is  commonly  carried  out  in  optics.  The 
expansion  (6)  in  terms  of  susceptibilities  applies  mutatis 
mutandis  to  the  transport  case,  and  to  our  knowledge  was 
hitherto  unknown. 

Depending  on  whether  F  is  polarized  in  the  x  or  z  direc¬ 
tion,  the  dominant  effect  will  differ;  for  F||x,  the  time  depen¬ 
dence  of  y(f)  will  dominate  through  the  ionization  rate, 
while  if  F||z,  e(t)  will  play  a  larger  role.19  Of  course,  e(f), 


I 

y(r),  and  T(t)  are  closely  interlinked  (see  comments  in  Ref. 
5),  but  for  our  purposes  we  can  treat  them  as  parameters 
determined  by  a  microscopic  theory  external  to  our  model. 

To  illustrate  the  different  types  of  behaviors  expected  in 
these  cases,  we  present  results  for  the  normalized  reflectivity 
(T'fr)2\R„(a  +  jLO\2  [(r/r,)2  =  0.017  is  the  peak  un¬ 
modulated  reflectivity]  of  a  high-quality  QW  in  Fig.  1.  We 
take  7i  =  0  and  consider  the  two  cases  2.5  (solid)  and 
1.52  ps-1  (dashed)  for  /x=0,  1,  and  2.  The  results  are  plot¬ 
ted  as  a  function  of  8=  Re  A.  [In  all  cases,  T  =  0.1  ps-1, 
y0=1.52  ps-1  (hy0-  1  meV),  and  £/27r=4  THz.]  This 
gives  the  reflected  probability  into  the  jxth  TS  as  a  function 
of  incident  frequency.  The  normalized  reflectivity  for  p=  1 
is  in  the  few  %  range  while  that  of  fx  —  2  in  the  10~4%  range 
for  the  parameters.  Figure  2  shows  Rtot(<*>)  for  the  same  pa¬ 
rameters.  The  relative  maximum  strength  of  the  first  TS  with 
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FIG.  3.  Total  reflected  probability  for  e^s/lir  ps-1  (dotted),  (l+i)ir  ps  1  (dashed),  and  isflir  ps  1  (solid).  The  parameters  are 
otherwise  those  of  Fig.  1 . 


respect  to  the  unmodulated  response  is  ~14.5%  for  ex 

-2.5  ps”1.  Figure  3  shows  Rtot(o))  for  ps”1 

(dotted),  (l+/)7r  ps”1  (dashed),  and  iyflTr  ps”1  (solid) 
with  the  other  parameters  the  same  as  in  Fig.  1  to  show  the 
effect  of  modulating  the  line  center  versus  the  linewidth.  In 

all  cases  | ps-1. 

We  turn  our  attention  to  time-domain  effects.  Specifically, 
we  describe  a  new  technique  for  coherent  control  of  excitons 
in  QW’s  using  phase-locked  pulse  pairs.  Previous 
schemes21,22  utilize  the  relative  phase  between  two  identical 
but  time-delayed  optical  pulses  coincident  with  a  narrow  op¬ 
tical  resonance  in  a  QW.  If  the  two  pulses  are  in  phase,  the 
interband  polarizations  excited  by  each  are  in  phase,  and  as  a 
result  ideally  four  times  the  population  of  excitons  associated 
with  a  single  pulse  are  generated;  if,  however,  the  pulses  are 
7 r  out  of  phase,  the  interference  is  destructive  and  the  QW  is 
coherently  depopulated  of  excitons  following  the  passage  of 
the  second  pulse.  Or  at  least  that  is  the  idea;  dephasing,  non- 
linearities,  and  inhomogeneous  broadening  all  conspire  to 
make  this  effect  less  marked  than  one  might  hope.22  In  the 
scheme  proposed  here,  two  identical  phase-locked  time- 
delayed  pulses  are  also  used,  but  rather  than  modulate  the 
relative  phase  of  the  two  pulses,  we  introduce  a  phase  shift  in 
the  polarization  induced  by  the  first  optical  pulse  by  a  half¬ 
cycle  THz  pulse  incident  on  the  QW  between  the  two  optical 
pulses.  Set  I(t)  =  I0(t  +  t0/2)+I0(t-t0f2)  where  70(O 
=  A(t)e~lEot\  A(t)  is  the  slowly  varying  envelope  of  a  single 
constituent  pulse  of  duration  r0<t0<yQ* ,  the  time  delay 
between  pulses.  Also  let  e(f)  =  e0  +  ei^(O  and  7(0“  7o- 
B(t)  is  a  peaked  function  of  time  centered  at  t-0  whose 
duration  is  much  less  than  t0.  For  example,  B(t)  might  be 
due  to  a  photoconductively  generated  half-cycle  THz  pulse 
propagating  down  a  transmission  structure  deposited  on  the 
QW.  For  t>t0 , 


T(t)  =  rf'°  dt'  e~^o(,~'')e~ie^i'd'"B(,")I(t') 

=  '*o'|>  ~ ie  i  f-  + 1  ]  |  dt'A(t')  (7) 

and  R(t)  =  T(t).  If  the  THz  pulse  is  chosen  so  that 
EtfZoodt" B(t")  =  (2n+  l)7r,  destructive  interference  be¬ 
tween  the  polarizations  excited  by  the  two  optical  pulses 
occurs  and  thus  coherent  depopulation  of  the  excitons  en¬ 
sues;  if  the  action  integral  is  2«7r,  then  constructive  interfer¬ 
ence  ensues. 

As  an  example,  if  the  half-cycle  THz  pulse  is  of  0.5  ps 
duration,  we  need  2  meV  to  achieve  a  i t  phase  shift. 
This  is  feasible  using  available  THz  techniques,  though  car¬ 
rying  out  this  experiment  may  be  difficult.  One  needs  very 
high-quality  QW’s  at  low  temperatures  to  ensure  narrow  ex- 
citonic  lines.  In  addition,  the  THz  transmission  structure 
should  be  designed  with  F(r)||z  to  suppress  modulation  of 
the  dephasing;  however,  this  can  be  compensated  for  by  de¬ 
creasing  the  amplitude  of  the  second  optical  pulse  in  the 
sequence.  Finally,  the  desirable  effect  will  be  degraded  if 
there  is  temporal  overlap  of  B(t)  and  A(t±t0/2),  thus  requir¬ 
ing  clean  optical  and  THz  wave  forms.  We  note  that  the 
present  scheme  in  no  way  circumvents  the  central  difficulties 
with  coherent  depopulation  of  the  exciton  level,  namely,  in¬ 
homogeneous  broadening  and  excitation-induced  dephasing. 

To  conclude,  we  present  an  exactly  solvable  model  for 
TSG  by  and  photon  propagation  through  THz-modulated 
QW’s.  Because  of  the  mapping  between  the  electromagnetic 
wave  equation  and  the  Schrodinger  equation,  light  propaga¬ 
tion  through  THz-modulated  QW’s  will  be  a  fruitful  way  to 
explore  resonant  tunneling  through  time-dependent  poten¬ 
tials.  We  believe  the  area  of  optical  analogues  of  time- 
dependent  quantum  mechanical  transport  problems — as  pio¬ 
neered  by  the  work  reviewed  in  Ref.  12 — is  in  its  infancy, 
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and  in  particular,  exploitation  of  optical  pulse  propagation 
through  QW’s  shows  great  promise  to  advance  the  field  mte- 
rially.  In  particular,  we  have  explored  the  effects  of  modu¬ 
lating  the  exciton  linewidth  and  energy  on  the  TS’s  appear¬ 
ing  in  the  optical  spectra.  We  have  also  proposed  an 


alternative  scheme  for  carrying  out  two-pulse  coherent  con¬ 
trol  of  excitons  in  QW’s. 

Helpful  discussions  with  T.  B.  Norris  are  gratefully  ac¬ 
knowledged.  This  work  was  supported  by  the  Office  of  Na¬ 
val  Research  and  the  Air  Force  Office  of  Scientific  Research. 
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clear  exciton  trapping  and  quasiadiabatic  population  transfer  in  a  two-band  semiconductor  is  predicted. 
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Driven  by  the  pursuit  to  understand  the  differences  and 
similarities  between  the  semiconductor  and  atomic  systems 
resonantly  excited  by  laser  pulses,  the  coherent  nonlinear 
dynamics  of  direct-gap  semiconductors  has  been  vigorously 
investigated  over  the  years.  Moreover,  with  the  ongoing  ad¬ 
vancement  of  short-pulse  laser  techniques  and  excellent 
semiconductors  samples,  classes  of  coherent  dynamic  phe¬ 
nomena  have  been  discovered  recently  including,  for  ex¬ 
ample,  coherent  exciton  control1,2  and  self-induced  trans¬ 
mission.3  Running  in  parallel  have  been  discoveries  in 
atomic  optics  such  as  population  trapping  using  frequency 
modulated  fields.  A  realization  of  such  state  trapping  and  its 
periodic  redistribution  in  a  two-level  atom  (TLA)  was  pre¬ 
dicted  in  Ref.  4.  In  the  past  year,  in  combination  with  other 
well-known  concepts  such  as  adiabatic  rapid  passage,  multi¬ 
photon  resonances,  and  Landau-Zener  transitions,  trapping 
was  experimentally  demonstrated  in  a  TLA  5.  For  two-band 
semiconductors,  however,  trapping  is  scarcely  expected  due 
to  Coulomb  many-body  complications  and  the  valence/ 
conduction  band-continua  of  free-carriers.  Surprisingly  in 
this  work  we  predict  clear  excitonic-state  trapping  (EST)  and 
quasiadiabatic  population  transfer  (QAPT)  using  multicol¬ 
ored  excitation  schemes. 

For  high-intensity  ultrashort  optical  pulses,  it  is  well  es¬ 
tablished  that  Rabi  oscillations  of  the  population  between 
two  states  can  be  seen  in  the  temporal  evolution  of  a  TLA.7 
However,  in  the  case  of  a  sinusoidal  frequency-modulated 
excitation,  square- wave  oscillations  of  the  population — 
periodic  state  trapping — as  well  as  more  complicated  and 
phase-dependent  structures  may  appear.  In  semiconductors, 
the  two-level  model  (as  a  first  approximation  to  a  two-band 
description)  is  considered  inappropriate  because  Coulomb 
many-body  interactions  result  in  a  renormalized  Rabi  energy 
and  bandedge,  and  excitation-induced  dephasing  (EID).  Nev¬ 
ertheless,  self-induced  transmission  and  multiple  Rabi  flop¬ 
ping  on  a  semiconductor  free  exciton  resonance  have  been 
recently  reported  in  bulk3  and  quantum  wells  (QW’s);6  these 
measurements  were  successfully  explained  within  the  frame¬ 
work  of  the  semiconductor  Bloch  equations.  Naturally  the 
question  arises:  Can  one  achieve  population  trapping  dynam¬ 
ics  in  semiconductors  using  frequency-modulated  or 
suitably-chirped  broadband  optical  pulses?  By  employing, 
for  example,  a  frequency-dependent  150-fs  full  width  at  half 
maximum  (FWHM)  irradiance  sech  pulse  excited  at  the  Is 
exciton  peak,  we  predict,  yes,  one  can  achieve  EST  and 
QAPT  in  a  semiconductor,  too.  We  present  trapping  dynam¬ 
ics  in  semiconductors  akin  to  trapping  in  atomic  systems.4 
We  discuss  in  detail  conditions  to  achieve  EST  and  QAPT 
even  when  a  broadband  pulse  excites  both  Coulomb-bound 
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electron-hole  (e-h)  pairs  (excitons)  and  free  e-h  pairs.  The 
other  well  known  trapping  phenomenon  in  atomic  systems  is 
the  coherent  population  trapping  effect,  which  has  been 
widely  used  inter  alia  in  adiabatic  population  transfer,  elec- 
tromagnetically  induced  transparency,  and  velocity  selective 
cooling.  As  pointed  out  recently  by  using  adiabatic  popula¬ 
tion  transfer  between  heavy-  and  light-hole  bands  (to  mimic 
a  three-level  atom),8  the  exact  analog  to  the  coherent- 
population-trapped  state  (dark  state)  used  in  atomic  physics 
is  not  possible  in  semiconductors. 

For  our  theoretical  approach  we  assume  the  validity  of  the 
rotating-wave  approximation  and  thus  neglect  the  possibility 
of  carrier-wave  Rabi  flopping.9  We  further  assume  a  two- 
band  QW  where  each  e-h  state  within  a  certain  band- 
structure  with  a  wave  number  k,  contributes  to  the  total  op¬ 
tical  polarization  P=2A~l'ZkdcvPkJ  with  dcv  the  transition 
dipole  matrix  element.  The  SBE  for  the  polarization  func¬ 
tions  are  (h=l)  (Ref.  10)  Pk=  WAkPk“/ftk(/k+/£--  1) 
+  Pk|corr  >  where  Ak  is  the  renormalized  energy  dispersion 
for  a  parabolic  two-band  semiconductor  and  flk  is  the  gen¬ 
eralized  Rabi  frequency,  while  f^h  are  the  carrier  distribu¬ 
tion  functions  of  the  electrons  and  holes.  The  Coulomb  in¬ 
teraction  is  treated  here  in  a  quasi-static  approximation 
(single  plasmon-pole).10  Similarly,  the  carrier  density  can  be 
calculated  from  the  electron  (or  hole)  population  N 
=  2A_12k/‘k,  where  /k  is  calculated  from  j^h=iP^i lk 
-*Pkflk +/kM|corr-  In  addition  to  the  terms  that  result  from 
the  time-dependent  Hartree-Fock  approximation,  the  carrier- 
carrier  (CC)  and  carrier-phonon  (CP)  collisions  drive  the 
nonequilibrium  distribution  functions  towards  quasiequilib¬ 
rium  Fermi  functions  and  yield  optical  dephasing.  The  CC 
collisions  (Coulomb  correlation  terms)  are  calculated  from 
the  e-h  Boltzmann  equations  including  also  a  correlation 
field,  nondiagonal  dephasing,  and  polarization  scat¬ 
tering.11-13  The  influence  of  CP  interactions  occurs  over 
much  longer  time  scales  than  those  studied  in  this  work  and 
can  be  safely  neglected. 

We  assume  input  optical  pulses  of  the  form  E(r ,t) 
—  E(t)e~ l[a)°t+k°z+  +  c.c.,  polarized  in  the  plane  of  the 

QW,  with  all  material  parameters  corresponding  close  to 
InxGaj _x As/GaAs  QW’s  (Ref.  14)  and  a)0-u)x  (Is  exciton 
resonance).  I^Ga^As/GaAs  QW’s  are  advantageous  since 
with  compressive  strain  one  can  increase  the  splitting  of  the 
heavy-  and  light-hole  exciton  and  thus  neglect  the  light-hole 
states,  therefore  validating  the  two-band  model.  In  analogy 
with  the  atomic  systems  we  modulate  the  input  field  by  the 
phase  factor  —  M  sin(ffr),  with  M  and  fi  the  index  and 
frequency  of  modulation.  For  a  TLA  the  resonant  part  of  the 
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FIG.  1.  Irradiance  spectra  for  a  frequency-modulated  and  un¬ 
modulated  150  fs  optical  pulse.  Energy  is  in  units  of  the  exciton 
binding  frequency. 

interaction  Hamiltonian  vanishes  if  M  is  chosen  such  that  the 
bessel  function  Jq(M)~0  and  one  is  left  with  the  nonreso¬ 
nant  rapidly  oscillating  terms.4  The  solid-state  community 
may  recall  that  a  somewhat  similar  criterion  has  been  derived 
to  realize  dynamic  localization  in  infinite  lattices  driven  by  a 
harmonic  time-dependent  electric  field.15,16 

In  the  TLA  it  is  further  known  that  the  two  bare  levels 
cross  at  the  times  tn  =  n  tt/20  (n  =  integer).  The  crossing  of 
energy  levels  occurs  quite  naturally  in  the  instantaneous 
frame  of  the  modulated  field.  The  population  distribution  at 
the  crossings  can  be  understood  semiquantitatively  by  inte¬ 
grating  numerically  the  time-dependent  Schrodinger  equa¬ 
tion  or  by  employing  Landau-Zener  theory.4  For  semicon¬ 
ductor  the  situation  is  much  more  complex,  whereby  the 
optical  field  provides  coupling  between  two  bands,  and  a 
modulated  near-resonant  field  leads  to  coupling  of  higher 
spectral  components  of  the  field  to  larger  k  states  in  the 
band,  and  the  field  components  below  the  resonance  provide 
a  nonresonant  coupling  to  the  excitonic  transition.  Moreover, 
the  Coulomb  interaction  yields  exciton  and  plasma  induced 
many-body  effects.  These  conditions  prohibit  us  from  having 
a  simple  set  of  criteria  for  trapping.  Thus  we  propose  here  a 
novel  spectral  distribution  for  die  field,  which  has  its  central 
frequency  detuned  far  below  the  excitonic  resonance,  and,  by 
strategic  sweeping  of  the  instantaneous  carrier  frequency,  dy¬ 
namical  population  trapping  can  be  realized.  As  in  the  TLA 
case  we  do  not  necessarily  have  to  choose  M  to  be  strictly  a 
zero  of  the  Bessel  function. 

For  the  modulation  frequency  we  choose  O  =  8  meV  un¬ 
less  stated  otherwise.  As  mentioned  above,  for  the  unmodu¬ 
lated  field,  E(t),  we  employ  a  150-fs  pulse  excited  at  the  Is 
exciton  peak.  In  Fig.  1,  we  depict  the  spectral  irradiance  of 
the  unmodulated  and  frequency-modulated  fields.  With 
modulation,  a  much  larger  bandwidth  can  be  obtained  with  a 
series  of  larger  peaks  extending  far  below  the  exciton  reso¬ 
nance;  the  injected  pulse  profile  is  sufficiently  broad  to 
couple  many  different  excitation  modes.  One  recognizes  that 
the  dominant  spectral  components  are  well  below  the  band 
edge.  By  increasing  the  index  of  modulation  (A/),  the  spec- 


FIG.  2.  (a)  Pulse-induced  carrier  density  (in  units  of  inverse 
Bohr  radius  squared)  as  a  function  of  time.  The  chain  line  corre¬ 
sponds  to  M=14.9  but  with  a  frequency  modulation  of  0/2.  (b) 
The  polarization  density  corresponding  to  the  dashed  case  in  Fig. 
2(a)  at  several  snapshots. 

tral  width  of  the  pulse  increases.  In  contrast,  without  any 
modulation  the  spectrum  of  the  excitation  pulse  is  signifi¬ 
cantly  narrower  and  peaked  at  the  central  frequency  as 
shown  in  Fig.  4(b)  (solid  line;  see  later).  Although  such  a 
modulation  is  difficult  to  obtain  experimentally,  one  only 
requires  about  one  approximate  oscillation  to  obtain  the  ap¬ 
propriate  energy  level  crossings.  Indeed  one  can  also  achieve 
trapping  phenomena  using  suitably  chosen  linear  chirps  (dis¬ 
cussed  below).  Further,  recent  developments  of  free-electron 
lasers  as  well  as  terahertz  solid-state  emitters  may  be  utilized 
in  materials  to  achieve  the  desired  frequency  modulation  in  a 
suitably  chosen  nonlinear  material.17 

For  an  optical  pulse  that  has  roughly  an  area  (integral  of 
its  Rabi  frequency  over  time)  of  67t  (~1  GW/cm2),  Fig. 
2(a)  shows  the  excitation-induced  density  for  the  QW  opti¬ 
cally  excited  at  the  Is  exciton  resonance  using  (i)  an  un¬ 
modulated  input  pulse  [the  solid  curve  depicts  the  carrier- 
density  showing  the  familiar  Rabi  flopping7,18  of  the  density 
(complete  inversion  is  not  possible  due  to  Coulomb  pro¬ 
cesses)];  and  (ii)  a  frequency-modulation  of  the  input  pulse 
for  two  different  values  of  the  index  of  modulation  M\  (the 
dashed  and  dotted  curves  are  markedly  different  from  the 
unmodulated  case  and  display  strong  EST  and  ultrafast 
QAPT).  The  values  of  M  were  chosen  in  accord  with  the 
original  work  in  atoms,  while  the  pulse  area  was  chosen  to 
be  large  enough  to  drive  the  population  into  trapping-it  is 
also  close  to  values  employed  experimentally.6  During  the 
time  interval  from  about  -160  to  160  fs,  the  population  re¬ 
mains  trapped  in  the  exciton  level,  eventually  jumping  out  in 
the  wake  of  the  exciton-continuum  quasi-adiabatic  crossing. 
For  the  case  of  0/2  and  M- 14.9  (chain  curve),  the  trapping 
efficiency  is  much  less  as  expected  and  QAPT  does  not  oc- 
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cur,  since  the  crossing  occurs  just  once  and  the  on-resonance 
field  components  are  not  sufficient  to  create  the  trappinglike 
feature;  however,  even  in  a  regime  where  the  pulse  irradi- 
ance  is  negligible  there  is  some  evidence  for  small  density 
changes  at  around  360  fs.  In  all  cases,  there  are  signatures  of 
phase  interference  due  to  coherent  carrier  evolution  along 
different  interfering  pathways  and  EID.  We  would  like  to 
point  out  that  to  obtain  such  good  Rabi  flopping  one  must 
include  EID  at  a  microscopic  level  that  includes  both  nondi¬ 
agonal  dephasing  and  polarization  scattering;11  these  contri¬ 
butions  from  CC  scattering  reduce  the  interband  optical  line- 
widths  of  the  higher  k  states  and  thus  limit  the  higher-energy 
continuum  occupations  in  comparison  to  the  pure  dephasing 
(essentially  a  relaxation-time  approximation)  treatment. 
Physically,  this  is  important  since  the  leading  edge  of  the 
pulse,  which  is  detuned  below  the  exciton  resonance,  pre¬ 
pares  the  system  for  the  excitation  of  real  population.  A  state 
and  energy-independent  dephasing  time  would  result  in  erro¬ 
neous  large  dephasing  of  die  initially  virtual  excitation  and 
therefore  suppress  the  QAPT  and  EST.  Polarization  scatter¬ 
ing  is  also  known  to  result  in  a  transfer  of  oscillator  strength 
from  the  continuum  to  the  exciton.13 

To  highlight  the  difference  between  the  Rabi  flopping  and 
trapping  we  depict  in  Fig.  2(b)  the  polarization  density 
(wave-packet),  |P(r,r)|2  with  r=re—rh  at  various  temporal 
snapshots.19  During  the  period  -260  to  100  fs  the  population 
becomes  strongly  trapped  in  the  exciton  state  indicated  in  the 
figure  by  a  high  probability  of  finding  the  electron  and  hole 
at  the  same  relative  position.  However  at  the  later  time  of 
280  fs,  after  the  density  increases  rapidly  (we  mention  that 
this  increase  is  almost  steplike  in  the  absence  of  dephasing), 
the  excitonic  probability  decreases  substantially  and  the 
wave  packet  speads  out  significantly,  demonstrating  that  the 
population  is  no  longer  trapped.  This  drammatic  spreading  of 
the  wave  packet  arises  due  to  the  modulation  and  the  result¬ 
ing  crossing  of  carriers  into  the  continuum.  The  above  pic¬ 
ture  sheds  much  more  light  on  the  trapping  scenario  than,  for 
example,  f^h  or  k,  which  do  not  clearly  distinguish  between 
Coulomb-bound  excitons  and  free  carriers  in  the  continuum. 
The  spatial  polarization  dynamics  at  r-0  may  in  fact  be 
probed  experimentally  using  conventional  four-wave  mixing 
techniques.  We  observe  here  a  method  to  control,  coherently, 
the  excitonic  wave  packet  by  tailoring  the  conditions  of  the 
energy  level  crossing.  To  probe  the  entire  spatial  dynamics 
one  would  need  to,  e.g.,  couple  a  THz  field  with  the  optical 
field;  this  would  involve  a  significantly  more  complex  analy¬ 
sis.  We  also  mention  that  to  direcly  detect  density  changes, 
experimentally  pulse-propagation  studies  are  very  difficult 
and  one  should,  for  example  ,use  the  technique  reported  re¬ 
cently  in  Ref.  6.  This  allows  for  the  detection  of  density 
changes  via  simple  differential  transmission  changes  of  a 
probe  pulse. 

The  aforesaid  phenomena  beg  a  more  physically  intuitive 
explanation:  Adding  a  frequency-modulation  to  the  excita¬ 
tion  field  is  equivalent  to  modulating  the  energy  separation 
between  the  ground  and  exciton/continuum  states.  In  Fig.  3 
we  show,  schematically,  the  various  crossings  of  energy  lev¬ 
els  that  lead  to  quasiadiabatic  transfer  of  population  at  these 
crossings.  The  crossings  “R-D”  transform  into  anti¬ 
crossings  (solid  lines)  due  to  the  coupling  with  the  effective 
Rabi  field.  In  contrast  the  crossing  at  “A”  is  unaffected  due 
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FIG.  3.  Schematic  of  the  energy  level  crossing  that  results  from 
modulation  of  the  field,  where  “A-D”  represent  the  various  cross¬ 
ings.  The  arrows  depict  the  evolution  of  the  population  (large  black 
dots). 

to  the  weak  coupling  of  the  v<->c  transition  on  the  leading 
edge  of  the  pulse.  The  modulation  of  the  dressed  energy 
levels  results  from  transforming  into  a  frame  corresponding 
to  the  instantaneous  field  frequency,  and  are  given  as  £c 
=  -[Aj;-A1],  2X=0  and  2V=  “[A1  +  cos(flr)],  for 
the  electrons  in  the  conduction  band  (c),  excitonic  state  (*), 
and  electrons  in  the  valence  band  ( v),  respectively.  One  can 
of  course  depict  the  diagram  in  reverse  for  the  holes.  Here, 
A1  (A£)  denote  the  detuning  of  the  central  laser  field  fre¬ 
quency  from  the  excitonic  (continuum)  energy  levels.  The 
coupling  of  the  field  to  the  v^x  and  v<->c  transitions  trans¬ 
forms  these  crossings  into  avoided  crossings.  Initially  the 
system  will  be  off  resonance  and  far  from  any  crossing;  as 
the  modulation  changes,  the  population  is  swept  though  reso¬ 
nance  and  it  evolves  quasi-adiabatically  into  the  excitonic 
state  where  it  displays  the  trapping  feature.  The  system  fur¬ 
ther  encounters  two  closely  spaced  crossings,  resulting  in  an 
enhanced  step-like  transfer  of  population  into  the  continuum. 
The  arrows  in  Fig.  3  indicate  one  possible  temporal  path 
(high  probability)  along  which  the  electrons  may  evolve  (AJ; 
is  fixed  for  simplicity).  The  detuning  term  Ajj-A1  is  the  dif¬ 
ference  between  the  continuum  energy  levels  and  frequency 
of  the  the  Is  exciton  peak.  One  should  keep  in  mind  how¬ 
ever  that  the  above  model  is  grossly  simplified  and  many- 
body  effects,  included  in  our  numerical  results,  will  compli¬ 
cate  things  substantially;  however  our  quantitative 
theoretical  study  is  in  fairly  good  agreement  with  the  above 
level  crossing  model. 

The  essential  difference  in  the  trapping  criterion  between 
the  atomic  case  and  semiconductor  case  lies  in  the  frequency 
content  of  the  exciting  field.  In  the  atomic  case  the  trapping¬ 
like  phenomenon  with  nearly  complete  inversion  results 
from  correlated  sideband  excitation  of  the  atom.20  The  fre¬ 
quency  components  of  the  modulated  field  excite  the  atom 
symmetrically  about  the  atomic  resonance,  resulting  in  trap¬ 
ping.  In  the  semiconductor,  a  symmetric  frequency  content 
of  the  excitation  does  not  lead  to  the  desired  trapping  due  to 
its  nonsymmetric  coupling  to  the  band  structure  of  the  semi¬ 
conductor;  effectively  the  high-frequency  components  of  the 
modulated  pulse  selectively  excite  the  continuum  of  states 
resulting  in  large  carrier  generation,  thus  washing  out  the 
trapping.  We  circumvent  this  problem  by  employing  a  modu¬ 
lated  pulse  with  its  predominant  frequency  content  away 
from  the  band  edge. 
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FIG.  4.  (a)  Pulse-induced  carrier  density  as  a  function  of.  (b) 
Corresponding  input  pulse  spectra  for  (a),  (c)  The  polarization  den¬ 
sity  for  the  positive  chirp  corresponding  to  the  dashed  case  in  (a)  at 
several  snapshots  (see  text),  (d)  As  in  (c)  but  for  negative  chirp. 


Next,  we  employ  linear-chirped  input  pulses  to  discern  if 
some  of  the  same  qualitative  trapping  features  (EST  and 
QAPT)  can  also  be  obtained  using  suitably  chosen  chirp  pa¬ 
rameters.  To  best  immitate  the  frequency-modulated  spec¬ 
trum,  we  choose  a  frequency  chirp  of  the  form:  <w— mux 


-120  m eV+acf  with  ac  being  ±0.5  meV/fs.  E0(t)  is  the 
same  as  before.  The  chirped-pulse  spectrum  along  with  the 
unchirped  spectrum  is  displayed  in  Fig.  4(b).  The  pulse- 
induced  density  with  positive  ac  is  shown  by  the  dashed  line 
which,  although  does  not  show  signs  of  trapping,  does  show 
fast  oscillations  in  the  density  as  well  as  a  rapid  increase  in 
the  density  near  the  crossing  time.  For  the  negative  ac, 
strong  excitonic  trapping  is  again  achieved  and  seems  to 
maintain  its  trapped  state  since  no  further  crossing  takes 
place.  This  is  again  clear  from  the  simplified  zero-order 
energy -level  picture  with  appropriate  crossing  of  the  1 5  ex- 
citon  state,  ground  crystal  state,  and  continuum  states.  The 
slope  and  sign  of  the  chirp  ( ac )  gives  us  a  handle  to  selec¬ 
tively  suppress  or  enhance  the  excitonic  transition.  The  trap¬ 
ping  and  population  transfer  are  highlighted  in  Figs.  4(c)- 
4(d)  which  show  the  wave-packets  at  the  times  -80  fs  (solid 
curve),  100  fs  (dashed  curve),  and  280  fs  (dotted  curve)  for 
both  the  negative  and  positive  chirp  pulse  excitation.  Popu¬ 
lation  transfer  is  not  expected  for  the  negative-chirp  case  as 
the  zero-order  energy  level  never  cross. 

In  conclusion,  we  predict  the  possibility  of  exciton  trap¬ 
ping  and  quasiadiabatic  population  transfer  in  a  two-band 
semiconductor  using  frequency-modulated  optical  pulses.  It 
is  clear  that  the  trapping  feature  arises  out  of  an  interplay  of 
the  dominant  excitonic  resonance  and  the  excitation  by  a 
sufficiently  strong,  broadband  off-resonant  pulse ,  with  a 
weak  yet  broad  spectral  content  to  excite  the  continuum  of 
states.  The  population  redistribution  results  from  the  crossing 
of  energy  levels.  Besides  being  an  intriguing  theoretical 
study,  our  results  are  timely  with  recent  advances  in 
frequency-modulated  spectroscopy  techniques  and  the  obser¬ 
vation  of  multiple  Rabi  flopping  on  free  exciton  transitions. 
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expanding  rurrner  tne  commercial 
applications  of  QPM  devices.2  How¬ 
ever,  QPM  engineering  finds  novel 
important  applications  far  beyond 
pure  frequency-conversion  devices. 

We  have  recently  exposed  the  great 
promise  of  engineered  QPM  patterns 
in  soliton  systems.  Solitons  in  QPM 
structures  form  by  the  mutual  para¬ 
metric  trapping  of  the  multi-fre¬ 
quency  light  waves  that  propagate  in 
the  crystal.3  The  “glue”  that  keeps  the 
multi-frequency  signals  locked  and 
trapped  in  robust  soliton  light  states 
is  their  energy-  and  phase-exchanges, 
mediated  by  cascading  of  quadratic 
nonlinearities.  Therefore,  all  the  sta¬ 
tionary  and  dynamic  features  of  the 
soliton  signals  depend  crucially  on 
the  local  phase-mismatch  experi¬ 
enced  by  the  multi -frequency  waves 
that  form  the  solitons. 

The  breakthrough  introduced  by 
QPM  engineering  is  the  possibility  to 
tailor  the  longitudinal  and  transverse 
variations  of  the  phase-mismatch 
between  the  waves,  hence  to  tailor  the 
properties  of  the  solitons.  This  is  a 


unique  possibility  that  links  the  implementation  of  poten¬ 
tial  soliton  devices  with  the  identification  and  design  of 
innovative  topologies  of  QPM  patterns.  In  particular,  we 
have  shown  the  soliton  based  signal  compression  and  gen¬ 
eral  shaping,  which  includes  not  only  tailoring  the  soliton 
width  and  profile  but  also  the  fraction  of  energy  carried  by 
each  of  the  waves  forming  the  multi-color  signals,  in  QPM 
structures  with  engineered  longitudinal  chirps;4  and  the 
spatial  soliton  switching  in  QPM  settings  with  engineered 
transverse  variations  of  the  phase  mismatch.5  Transverse 
patterns  illustrative  of  such  possibilities  include  chirps,  tilts, 
dislocations  and  wells,  as  those  shown  in  Figure  la.  In  the 
particular  example  shown  in  Figure  lb,  upon  strong  illumi¬ 
nation  the  soliton  crosses  the  QPM  well.  However,  with  a 
weaker  illumination  the  soliton  is  trapped  in  the  well.  Soli¬ 
ton  self-bouncing  and  power-dependent  soliton  scanning 
are  other  potential  feasible  possibilities  shown.  More 
advanced  and  practical  QPM  lay-outs,  including  modulat¬ 
ed,  multi-periodic,  and  quasi-periodic  gratings,  and  engi¬ 
neered  QPM  multi-grating  arrays  are  now  under 
investigation.  Engineered  competing  quadratic  and  cubic 
effective  nonlinearities,6  bandwidth-tailored  soliton  excita¬ 
tion,  quasi-discrete  soliton  signals,  and  novel  kinds  of  stable 
solitons  constitute  important  examples  of  the  fascinating 
potential  outcome. 

We  believe  that  these  findings  constitute  the  tip  of  the 
iceberg  of  the  potential  of  QPM  engineering  for  quadrat¬ 
ic  soliton  control  devices  in  single-,  multi-pass,  and  cavity 
geometries.  We  anticipate  that  creative  QPM  engineering 
might  offer  crucial  building  blocks  of  all  future  photonic 
devices  employing  quadratic  soliton  signals. 
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Terahertz  Sources 


Extremely-Tunable  Terahertz  Emission:  Coherent 
Population  Hopping  in  a  DC-Biased  Quantum  Well 

S.  Hughes  and  D.S.  Citrin 

Terahertz  (THz)  light-emitters  are  a  welcome  addition 
to  the  laser  engineer  s  palette  since  they  offer  a  marriage 
of  optics  and  high-speed  electronics.  In  the  field  of  semi¬ 
conductor  nanostructures,  the  exploration  of  THz  genera¬ 
tion  schemes  and  THz- induced  carrier  dynamics  provides  a 
useful  means  to  study  charge-carrier  transport.  The  inter- 
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play  of  optical  and  THz 
transients,  moreover, 
provides  a  unique  win¬ 
dow  for  simultaneously 
studying  inter-  and 
intra-band  carrier 
dynamics  in  semicon¬ 
ductor  heterostructures. 
Applications  with  the 
semiconductor  laser1 
include  gain  modulation 
and  all-optical  ultrafast 
switching  rates  of 
around  500  gigabit/ 
sec — a  useful  prerequi¬ 
site  for  the  next  genera¬ 
tion  of  broad-band¬ 
width  optical  communi¬ 
cations  networks  and 
the  inevitable  gigabit 
Internet.  More  funda¬ 
mental  applications 
encompass  coherent 
control  of  electronic 
wavepackets,2  nanotech¬ 
nology,  and  quantum 
computing.  In  short  the 
cry  for  reliable  THz 
emitters  can  be  heard 
throughout  the  globe. 
However,  apart  from 
expensive  free-electron 
lasers,  a  common  limita¬ 
tion  to  the  generation  of 
THz  transients  in  the 
laboratory  is  the  lack  of 
tunable  sources. 

The  successful  advancement  in  short-pulse  laser  tech¬ 
niques  and  excellent  semiconductors  samples  led  to  a  break¬ 
through  in  the  last  year  by  the  clear  observation  of 
self-induced  transmission  (SIT)  in  bulk  semiconductors3 
and  multiple  Rabi  floppping  in  quantum  wells  (QWs).4  Rabi 
flopping  is  a  well  known  coherent  effect  in  atomic  systems 
but  has  long  been  considered  to  be  the  Holy  Grail  to  the 
ultrafast  semiconductor  community — who  are  more  than 
ever  driven  by  the  pursuit  to  understand  the  differences  and 
similarities  between  the  semiconductor  and  atomic  systems 
resonantly  excited  by  laser  pulses.  Very  recently  we  have 
devised  a  new  method5  that  exploits  Rabi-flopping-induced 
density  oscillations  in  a  weakly-DC-biased  QW  to  realize 
very  intense  (~kVs/cm),  tunable  terahertz  transients.  In  our 
scheme,  extreme —  tailored — tunability  is  achieved  by  simply 
varying  the  area  of  the  input  optical  pulse.  We  employ  a  150 
fs  FWHM  ir radiance  pulse.  For  a  weakly-biased  QW  we 
assume  d  =  0.8  ^(corresponding  bias  field  ~  10-100  Vic m), 
and  calculate  the  emitted  field  at  r  =  1  cm.  We  assume  the 
validity  of  the  slowly  varying  envelope  approximation 
(which  is  not  always  valid,  see  Reference  6),  and  compute  the 
semiconductor  Bloch  equations  including  all  the  relevant 
scattering  mechanisms. 

Figure  la  shows  the  excitation-induced  density  for  the 


QW  optically  excited  at  the  Is  exciton  resonance  using 
the  input  Rabi  energies  of  15  meV  and  30  meV  (pulse 
irradiance:  ~  GW/cm2),  respectively.  The  curves  display 
the  conduction  band-valence  band  Rabi  flopping3,4 
(complete  inversion  is  not  possible  due  to  Coulomb 
processes).  Material  parameters  correspond  to 
I^Ga^As/GaAs  QWs  and  agree  well  with  the  experi¬ 
mental  observations.  In  Figure  lb  we  show  the  emitted 
THz  field  versus  time  corresponding  to  the  density  oscil¬ 
lations  in  Figure  la.  Each  transient  is  approximately  500 
fs  in  duration  (due  to  dephasing)  and  undergoes  a  series 
of  polarity  changes.  In  the  EM  spectra  (Figure  lc),  clear 
peaks  in  the  THz  regime  are  seen.  For  larger  areas,  more 
Rabi  flops  occur3,5  that  subsequenly  lead  to  THz  emis¬ 
sion  with  larger  frequencies.  We  note  that  current  tech¬ 
niques  employing  sub-ps  optical  pulses  for 
THz-generation  are  usually  limited  to  about  the  1-3 
THz  regime  at  best.  Our  selected  examples  here  easily  go 
to  the  lOs-of-THz  regime  with  the  added  bonus  of  much 
larger  fields.  Futhermore  they  do  not  suffer  from  any 
bandwidth  limitations  (c.f.  optical  rectification  tech¬ 
niques).  These  studies  are  motivated  in  part  because  of 
recent  dramatic  advances  in  high-speed,  hybrid  light¬ 
wave -THz  photonic  applications. 
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Real-Time  Optical  Pulse  Characterization  Using  SPIDER 

Matthew  E.  Anderson,  Timothy  M.  Shuman,  Chris  laconis  and 
Ian  A.  Walmsley 

The  recent  advances  in  ultrafast  lasers  have  mandated  a 
change  in  the  techniques  used  to  measure  them.  It  is 
no  longer  sufficient  to  measure  merely  the  spectrum  or 
the  auto-correlation  of  an  ultrafast  laser  pulse;  a  more 
complete  picture  is  required,  and  this  involves  determin¬ 
ing  the  phase  of  the  time-varying  optical  field.  The  ability 
to  rapidly  change  the  shape  of  the  pulse  is  a  key  technolo- 
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Hughes  Figure  1.  (a)  Temporal  variation  of 
carrier  pulse-induced  carrier  density  in  units  of 
inverse  Bohr  radius  squared;  excitation  is  at 
the  Is  exciton  resonance.  The  peak  amplitude 
of  the  pulse  Is  15  meV  (solid  curve)  and  30 
meV  (dashed  curve) .  (b)  Emitted  THzfield  ver¬ 
sus  time  for  the  15  meV  and  30  meV  pulses  at 
r*  1  cm  from  the  excitation  spot,  (c)  Irradi¬ 
ance  of  the  emitted  THz  spectra  corresponding 
to  (b). 
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